Entanglement and Quantum Phase Transition in Low Dimensional Spin Systems 



in 
o 
o 

<N 



Yan Chen\ Paolo Zanardi^, Z. D. Wang^-^, and F. C. Zhang^ 

^Department of Physics and Center of Theoretical and Computational Physics, 
University of Hong Kong, Pokfulam Road, Hong Kong, China 
'^Institute for Scientific Interchange Foundation, Viale Settimio Severo 65, 1-10133 Torino, Italy 
^National Laboratory of Solid State Microstructures, Nanjing University, Nanjing, China 

(Dated: February 1, 2008) 

Entanglement of the ground states in XX Z and dimerized Heisenberg spin chains as well as in 
a two-leg spin ladder is analyzed by using the spin-spin concurrence and the entanglement entropy 
between a selected sublattice of spins and the rest of the system. In particular, we reveal that 
quantum phase transition points/boundaries may be identified based on the analysis on the local 
extreme of this entanglement entropy, which is illustrated to be superior over the concurrence 
scenario and may enable us to explore quantum phase transitions in many other systems including 
higher dimensional ones. 
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Introduction. Quantum entanglement has attracted a 
lot of attention for its potential applications in quantum 
information processing 1, 2]. More recently, entangle- 
ment has also been recognized to play an important role 
in the study of quantum many particle physics, and ex- 
perimental measurements have demonstrated that it can 
affect the macroscopic properties of solids 3j . There have 
been a number of theoretical studies of the entanglement 
and quantum phase transitions Q in one dimensional 
spin systems and in interacting fermion and boson sys- 
tems UU,U,M^U,M,M,^M^MM,E^ These studies 
showed that the entanglement can exhibit a non-trivial 
behavior in condensed matter systems, such as the coin- 
cidence of the singularity in the entanglement and quan- 
tum phase transition point in certain systems 17J . Most 
of the previous works were mainly focused on spin chains 
and the spin-spin concurrence 1] was used to describe 
(two-particle) the entanglement. The scaling behavior of 
entanglement between a block of contiguous spins and 
the rest of the system in a spin chain has been studied 
both near and at the quantum critical point 18], and 
a connection with quantum phase transitions was elabo- 
rated. 

In this paper, we adopt a preferable and distinct way 
to partition the system to study the ground state of 
three low dimensional spin systems: the XXZ spin chain, 
dimerized Heisenberg chain, and two- leg spin ladder. We 
find that the entanglement entropy is superior over con- 
currence in revealing quantum transition points. In the 
dimerized Heisenberg system, this entanglement entropy 
has maxima and minima which are in one-to-one corre- 
spondence with the transition points, while the concur- 
rence fails to locate them. Even though the spin ladder 
system exhibits a complex phase structure, our entan- 
glement entropy scenario is likely to enable us to identify 
the phase boundaries. It is expected that our approach is 
more efficient and powerful in exploring quantum phase 
transitions in various systems. 

The models. We here study the entanglement in three 
types of spin-1/2 systems: XXZ chain, dimerized Heisen- 



berg chain, and two-leg XXZ ladder, with the Hamilto- 
nian given by 
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+ J' ^ {si^si^ + sy^^sy^^ + asi^si^) , (3) 

where Si is the spin operator, N is the linear dimension 
of the lattice (we adopt periodic boundary conditions). 
The anisotropy parameter is denoted by A in the XXZ 
model. Ji and J2 are two nearest-neighbor exchange cou- 
plings in the dimerized Heisenberg Hamiltonian. In the 
ladder system, the Hamiltonian contains intersite inter- 
actions along the chains (J) and the rungs (J'). 

To quantify the entanglement in these spin systems, 
we first adopt concurrence to describe two-spin entangle- 
ment Q. This quantity can be calculated from correla- 
tion functions G"" :=< a" (g) a" >, {a = x,y,z) |^: 
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where the two site indexes are omitted. Since the correla- 
tion function decays rapidly, the concurrence is nonzero 
only between two closer sites. Therefore, the concur- 
rence seems to provide limited information on the quan- 
tum phase transition. In a different approach, Vidal et 
al. 18] studied the entanglement between a block of L 
contiguous spins Bl and the rest of the system Rl- The 
entropy of entanglement for the ground state is given 

by, 



Sl ■■= -tr (pl log2 Pl) , 
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FIG. 1: The ground state concurrence (a) and the entangle- 
ment entropy per site Sl/L (b) in terms of A in XXZ spin 
chain for various lattice sizes [N — 2L). The inset of (b) 
displays the sublattice bipartition. 



where pL = trBL|^g)(^g| is the reduced density matrix 
for Bl^ a block of L spins. The scaling of the entan- 
glement in the block size L may allow one to establish 
a connection between Sl and the quantum phase tran- 
sition: S'l ~ InL at the critical point, and Sl ^ const 
away from the critical point; while this idea is applicable 
to few one dimensional models. Here we adopt a distinct 
way that bi-partitions the system into two sublattices Bl 
and Rl, with the largest number of connecting bonds be- 
tween them 0, 123 . The motivation of such partition is 
to best reveal the correlation between the two sublattices. 
We will see that in most cases the considered entangle- 
ment entropy between the two subsystems scales linearly 
with the size of the boundary between them [notice that 
in the case 18.1, t he boundary is zero-dimensional i.e., 
one single SDin]]21j. 

All the above three Hamiltonians have a rotational 
symmetry around the z-axes. The calculations presented 
below are carried out in an invariant subspace with 
5'^ = 0. Lanczos algorithms are employed to calculate 
the ground state |^^), from which we construct the den- 
sity matrix for the whole system. We then obtain the 
reduced density matrix Pl by tracing out Rl^ and com- 
pute its Von Neumann entropy in Eq.(5). 

XX Z -chain. Let us first consider the one-dimensional 
XXZ model. As a simple toy model, a great deal of work 
has been devoted to analyze its entanglement and quan- 
tum phase transition 22]. It is well known that A = 1 
is the antiferromagnetic isotropic point, A = — 1 is the 



ferromagnetic isotropic point, and A = is the pure XY 
point. This allows us to describe the various domains 
as a function of A. The system is in an Ising ferromag- 
netic phase at A < — 1, Ising antiferromagnetic phase at 
A > 1, and XY phase at -1 < A < 1. In Fig. 1(a), we 
show the pairwise concurrence as a function of anisotropy 
A for different number of sites L. Because the concur- 
rence here is expressed in terms of the two-site correlation 
function of the nearest-neighbor sites, its value converges 
quickly as L increases. As we can see, the concurrence at 
one phase transition point A = 1 reaches the maximum 
while the concurrence emerges at another transition point 
A = — 1 as L approaches infinity. These two quantum 
phase transition points can be identified from the analy- 
sis of concurrence in this case. We now address the en- 
tanglement entropy of the sublattice Bl = {odd — sites} 
with L = N/2 [23. Fig. 1(b) displays the numerical 
results of the entanglement entropy Sl/L as a function 
of A. We obtain the reduced subsystem of crosses by 
tracing out the spin degree of freedom on circle points. 
As the system is in the vicinity of the quantum phase 
transition point, we may expect Sl/L to reach its ex- 
treme value. We find that the transition points A = 1 
and A = — 1 correspond to the maximum and minimum 
of the sublattice entanglement entropy. As for the latter 
case, a simple analysis of the scaling behavior around the 
minimum point shows that the location of the transition 
point approaches to A = — 1 and Sl/L = as the size 
of the subsystem increases. Thus at L ^ cxo, Sl/L ^ 
for A < —1. Therefore, the two transition points can be 
clearly specified and a distinct connection between quan- 
tum phase transition and the entanglement entropy has 
been established. 

Dimerized Heisenberg chain. We now consider the 
dimerized Heisenberg chain |2j|; this model is charac- 
terized by an alternation of strong and weak bonds be- 
tween two nearest neighboring spins. In the case of 
< J2/J1 <C 1, the ground state is just an ensemble 
of N/2 uncoupled dimers around the strong bonds. Con- 
sequently, there is an energy gap of order of Ji to sep- 
arate the singlet ground state {Sz = 0) from the first 
excited state with Sz = ^- All the spins are locked into 
singlet states. At J2/J1 = 1, the system is reduced to 
the isotropic antiferromagnetic Heisenberg chain that is 
quasi-long-range ordered, and belongs to a different uni- 
versality class from the dimerized system. In the case of 
J2/ Ji < 0, the ground state is a product of spin singlet 
pairs coupled by ferromagnetic bonds. We can define 
two concurrences Ci and C2, which correspond to two 
nearest-neighboring spins coupled by bonds Ji and J2, 
respectively. The pairwise concurrence as a function of 
J2/J1 for different lattice size L is shown in Fig. 2(a). It 
can be easily seen that the value of concurrence is inde- 
pendent of the lattice size. In the limit of strong dimer- 
ization, J2/J1 = (or equivalently J2/J1 00), the 
concurrence Ci (C2) reaches a maximum value. However, 
the concurrence analysis does not enable us to identify 
the transition point J2/J1 = 1 unambiguously. Below, 




FIG. 2: For various lattice sizes A/" in a dimerized Heisenberg 
chain, (a) the ground state concurrences Ci (sohd hne) and 
C2 (dotted hne) versus A, and (b) Sl/L as a function of the 
sublattice size L. The inset in (b) shows the divergence of 
the dp {S L / L) / d{J2 / JiY ^he block entanglement minimum 
point as the system size increases. 



we examine the entanglement entropy of the chosen odd- 
site sublattice (or the enen-site sublattice). In this case, 
the two transition points are found to correspond to a 
local maximum and a local minimum in the entangle- 
ment entropy [see Fig. 2(b)]. Here it is interesting to 
notice that the second derivative of Sl/L with respect 
to J2/J1 seems to exhibit singularity at the minimum 
point J2/J1 = 1 (see Fig. 2(b) inset), which corresponds 
to the only gapless point. Due to the limitation of our 
computational resources, a detailed scaling analysis is not 
available. In addition, it is seen that Sl/L is weakly size 
dependent and converges quickly to a fixed value at the 
transition point with increasing lattice size. 

Spin ladders. Spin ladders have been attracting great 
interest in the field of low dimensional spin systems in 
recent several years The physics of Heisenberg 

ladders can be understood in some special limits. For 
large and negative J', the 2-leg ladder is equivalent to 
an 5* = 1 spin chain, and the system is gapful 26]. For 
large and positive J', the two spins along the rung are 
locked into a singlet state. The ground state is a collec- 
tion of spin singlets on each rung and is gapful. The phase 
diagram of two coupled XXZ chains was previously de- 
termined by utilizing Hamiltonian mappings and Abelian 
bosonization 27]. Previous studies presumably indicated 
five distinct phases in this complicated phase diagram as 
a function of A and J' /J. We first examine the depen- 
dence on the coupling constant of two nearest-neighbor 
spin concurrence: this analysis failed to describe the rich 



structure of the phase diagram. Below we report our 
study on the entanglement entropy of the chosen sub- 
system and its relevance to quantum phase transition 
point. The spatial profiles and contour plots of Sl/L 
as a function of J' / J and A are displayed in Fig. 3. 
From our analysis from above-mentioned two models, we 
know that both the ridges and valleys may correspond to 
possible phase boundaries. For example, one can distin- 
guish the boundary J'/J = (ridge or valley), A = 1 for 
jyj < (ridge), A = -1 for J' < (ridge). Derived 
from the numerical results. Fig. 4 shows the schematic 
phase boundaries of a coupled 2-leg XXZ spin ladder sys- 
tem (the dashed line denotes the ridge of its derivative). 
Most phase boundaries are in agreement with those in 
Ref. 26. Our results suggest the existence of a new quan- 
tum phase in the region A < —1.5 and J' / J > 0.5, which 
has not been reported in previous studies [23, HE • 

Finally, it is important to emphasize that, in compari- 
son with the previous investigations 118], we do not have 
to deal with a large system i.e., a large L, to capture its 
scaling behavior. This indicates that the finite size effects 
are not important in our approach; numerical results for 
a quite small system may disclose relevant information 
about the quantum phase transition points. In this sense, 
the present approach is broadly applicable to many other 
systems including spin and electron systems in higher di- 
mension 19], not limited merely to the spin-chains. 

Summary. We have investigated the ground-state en- 
tanglement in three one-dimensional spin systems by 
using both two-spin concurrence and entanglement en- 
tropy of the preferably chosen sublattice. Concurrence 
may provide some partial insights about certain quan- 
tum phase transition points. However, our entanglement 
entropy scenario allows us to establish a distinct con- 
nection between its local maxima/minima and transition 
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points, which is promising for shedding a new hght on the 
understanding of quantum phase transition and quantum 
entanglement. 
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